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In  an  opposed-jet  diffusion  flame  experiment,  under  certain  conditions,  after  the  extinction  of  the  dif¬ 
fusion  flame,  an  edge  flame  can  be  obtained.  This  ring-shaped  edge  flame  was  first  reported  in  1959  by 
Potter  and  Butler  but  received  little  attention.  It  was  reported  again  recently  in  a  numerical  and  an 
experimental  work  and  is  responsible  for  an  interesting  transition  between  two  distinct  burning  flames 
(multiple  solutions).  Motivated  by  our  previous  numerical  results,  obtained  with  simplified  kinetics  and 
some  recently  reported  experimental  data,  we  performed  direct  numerical  simulations  of  this  transition  to 
investigate  the  underlying  physical  mechanisms.  The  appearance  of  an  edge  flame  after  the  extinction  of 
the  diffusion  flame,  the  hysteresis  reported  in  the  experiments,  and  the  existence  of  multiple  vigorously 
burning  flames  at  identical  conditions  are  all  captured  by  our  simulations.  Our  numerical  results  show  that, 
in  the  absence  of  an  inert  coflow  curtain,  when  the  diffusion  flame  disk  is  extinguished,  an  edge  flame 
forms  and  propagates  in  the  mixing  layer.  After  the  formation  of  this  edge  flame,  even  when  the  applied 
strain  rate  is  reduced  to  the  initial  suberitical  value,  the  diffusion  flame  disk  does  not  reappear,  because 
the  local  fluid  velocity  still  exceeds  the  propagation  speed  of  the  edge  flame.  This  hysteresis  has  significant 
implications  in  the  common  submodel  that  utilizes  the  strain  rate  as  a  parameter  to  determine  local 
reignition  in  flamelet  models;  it  indicates  that  a  suberitical  strain  rate  is  not  a  sufficient  condition  for  the 
reignition  of  a  diffusion  flame.  Further  investigation  of  this  phenomenon  is  clearly  needed  to  refine  sub¬ 
models  of  local  extinction  and  reignition  in  the  flamelet  models  for  turbulent  diffusion  flames.  The  opposed- 
jet  configuration  provides  a  convenient  platform  to  analyze  edge  flames  which  are  stabilized  aerodynami- 
cally  in  a  two-dimensional  geometry,  thus  making  matching  two-dimensional  direct  numerical  simulations 
effective. 


Introduction 

We  report  direct  numerical  simulations  (DNSs)  of 
a  transition  of  flame  configuration  recently  reported 
in  the  experiments  of  Pellett  et  al.  [1]  on  a  diffusion 
flame  of  diluted  hydrogen  and  air  stabilized  on  an 
opposed-jet  burner.  A  similar  phenomenon  was  re¬ 
ported  in  one  of  our  early  papers  [2],  where  we  an¬ 
alyzed  the  effects  of  the  Damkohler  number  on  the 
flame  configuration  in  an  opposed-jet  geometry.  In 
fact,  as  early  as  1959,  in  the  pioneering  work  of  Pot¬ 
ter  and  Butler  [3],  such  a  ring-shaped  edge  flame 
was  observed;  it  was  then  called  the  “broken”  flame. 
Even  though  the  opposed-jet  configuration  is  used 
extensively  in  studying  strained  diffusion  flames,  the 
transition  involving  an  edge  flame  is  rarely  reported, 
primarily  because  such  edge  flames  exist  after  the 
local  extinction  of  the  strained  diffusion  flame  only 
when  there  is  no  inert  gas  coflow  curtain — a  very 


common  practice  in  experiments  to  reduce  the  ef¬ 
fects  of  air  entrainment  and  buoyancy  and  to  elimi¬ 
nate  the  “edge”  of  the  flame  in  order  to  create  a  disk¬ 
shaped  flame  and  facilitate  laser  diagnostics. 

Tomboulides  et  al.  [2]  reported  a  two-dimensional 
simulation  of  an  opposed-jet  diffusion  flame  using 
one-step  chemistry  and  simplified  transport  prop¬ 
erties  without  the  inert  coflow  curtain.  It  was  found 
that  when  the  Damkohler  number  was  increased  be¬ 
yond  a  critical  value,  the  flame  structure  changed 
from  a  disk-shaped  diffusion  flame  to  an  edge  flame. 
It  was  shown  that  the  edge  flame  was  a  tribrachial, 
or  triple,  flame  with  rich  and  lean  premixed  arms 
and  a  trailing  diffusion  flame.  In  Ref.  [1],  the  large 
hysteresis  between  “blowoff  ’  and  “restoration”  (ter¬ 
minology  of  Ref.  [1])  and,  in  particular,  the  descrip¬ 
tion  of  the  ring-shaped  flame  resemble  what  was  re¬ 
ported  in  Ref.  [2],  We  suspect  the  underlying 
physical  mechanism  for  the  phenomena  observed  in 
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Ref.  [1]  is  the  formation  and  propagation  (with  posi¬ 
tive  propagation  velocity)  of  an  edge  flame.  Unfor¬ 
tunately,  no  further  information  on  the  structure  of 
the  edge  flame,  such  as  flame  shape  visualization, 
temperature,  and  so  forth,  are  available  in  Ref.  [1], 
The  similarity  between  the  experimental  findings  [1] 
and  our  earlier  numerical  results  [2],  the  novelty  of 
using  an  opposed-jet  setup  to  study  edge  flames,  and 
its  implications  on  the  modeling  of  extinction  and 
reignition  of  a  diffusion  flamelet  [4,5]  prompted  us 
to  investigate  this  diffusion/edge  flame  transition 
further  using  detailed  two-dimensional  transient 
simulations  with  conditions  that  match  those  used  in 
the  experiment. 

In  general,  an  edge  flame  exists  on  the  boundary 
of  a  diffusion  flame.  Depending  on  the  sign  of  its 
propagation  velocity,  it  can  be  an  ignition  wave  or  a 
failure  wave  [6],  However,  unlike  a  premixed  flame, 
the  propagation  speed  of  an  edge  flame  is  dependent 
on  many  parameters  in  addition  to  the  thermody¬ 
namic  properties  (including  kinetics  and  transport 
properties)  of  the  mixture.  Its  ability  to  propagate 
with  a  positive  or  negative  velocity  is  responsible  for 
the  spread  or  blowout  of  a  diffusion  flame.  It  also 
offers  the  mechanism  for  a  diffusion  flame  to  anchor 
itself  [6-11],  Since  an  edge  flame  exists  whenever  a 
diffusion  flame  has  a  boundary  (such  as  a  hole  on  a 
diffusion  flame  sheet),  it  is  known  to  play  an  impor¬ 
tant  role  in  turbulent  diffusion  flames,  as  it  is  re¬ 
sponsible  for  the  local  phenomena  of  extinction  (fail¬ 
ure  wave)  and  reignition  (ignition  wave)  [4,10,12], 
The  triple  flame,  a  special  case  of  edge  flames,  has 
received  considerable  attention.  For  example,  the 
theoretical  work  of  Dold  [13]  and  Hartley  and  Dold 
[14],  the  one-dimensional  model  of  Buckmaster  [6], 
the  experimental/numerical  work  of  Kioniet  al.  [15], 
and  the  numerical  work  of  Ruetsch  et  al.  [16]  all 
illustrate  that  the  propagation  velocity  of  the  triple 
flame  depends  on  the  transverse  gradient  of  the  mix¬ 
ture  concentration,  flame  curvature,  and  transport  in 
the  transverse  direction.  Two-dimensional  numerical 
simulations  of  triple  flames  in  the  coflow  geometry 
for  both  the  transient  and  steady  (burner  stabilized) 
cases  with  chemical  kinetics  and  transport  models  of 
various  complexities  have  recently  been  reported 
[17-20], 

In  this  paper,  we  focus  on  the  multiple  vigorously 
burning  (as  opposed  to  the  weakly  burning  or  pure 
mixing  solution)  stable  flames  observed  in  an  op¬ 
posed-jet  configuration  [1],  We  perform  calculations 
with  conditions  that  match  one  set  of  experimental 
conditions  to  investigate  the  nature  of  the  two  dif¬ 
ferent  flames.  In  particular,  we  are  interested  in 
showing  that  an  edge  flame  is  indeed  responsible  for 
what  was  observed  in  Ref.  [1]  and  in  providing  ad¬ 
ditional  information  on  the  structure  of  this  edge 
flame.  Note  that  the  multiple  solution  phenomenon 
considered  here  is  unlike  the  familiar  S-shaped 
curve  with  the  upper  branch  representing  a  burning 


flame  and  the  lower  branch  corresponding  to  pure 
mixing  (or  a  weakly  burning  solution)  [6],  We  ob¬ 
served  two  stable  flames  in  addition  to  the  pure  mix¬ 
ing  solution — that  is,  three  stable  solutions.  To  the 
best  of  our  knowledge,  this  is  the  first  attempt  to 
perform  DNSs  of  multiple  vigorously  burning  solu¬ 
tion  branches  and  the  associated  hysteresis  effects. 
Finally,  we  would  like  to  note  that  in  using  the  op¬ 
posed-jet  geometry,  we  were  able  to  simulate  three 
different  types  of  edge  flames:  a  propagating  ignition 
wave  (positive  propagation  velocity  with  respect  to 
the  unburned  stream),  a  propagating  failure  wave 
(negative  propagation  velocity  with  respect  to  the 
unburned  stream),  and  also  the  flame  structure  at 
the  edge  of  the  diffusion  flame  disk.  Due  to  the 
length  limitation  on  the  symposium  papers,  these 
other  edge  flames  will  be  reported  elsewhere  in  the 
near  future. 


Methodology  and  Problem  Definition 

We  employed  a  direct  numerical  simulation  ap¬ 
proach  to  solve  the  transient,  two-dimensional,  low 
Mach  number  reduced  form  of  the  conservation 
equations  of  mass,  momentum,  energy,  and  species 
in  the  axisymmetric  cylindrical  coordinate  system. 
The  numerical  code  employed  is  based  on  the  spec¬ 
tral  element  method  for  spatial  discretization  and  a 
high-order  stiffly  stable  splitting  scheme  that  sepa¬ 
rates  the  thermochemistry  from  the  hydrodynamic 
subsystems  for  time  integration.  A  special  method 
was  employed  to  ensure  exponential  convergence  on 
the  axis  of  symmetry  [21],  Details  on  the  mathe¬ 
matical  formulation  and  the  numerical  method  can 
be  found  in  Refs.  [2,22],  while  some  of  our  previous 
studies  with  this  numerical  code  can  be  found  in 
Refs.  [23,24],  Recently,  we  implemented  the  matrix- 
free  stiff  ordinary  differential  equation  solver 
VODEPK  [25],  which  solves  the  thermochemistry 
part  much  more  efficiently  than  the  previous  solver 
(LSODES  [26]). 

Detailed  chemical  kinetics  [27]  and  realistic  trans¬ 
port  (the  mixture  fraction  formulation  with  a  correc¬ 
tion  diffusion  velocity  term  to  preserve  the  conser¬ 
vation  of  mass  [28])  were  utilized  in  the  simulations. 
Grid  independence  was  ensured  by  repeating  cal¬ 
culations  with  increasing  orders  of  interpolants,  and 
the  minor  effect  of  the  outflow  boundaries  on  the 
solution  was  verified  by  placing  them  at  increasing 
distance  from  the  axis  of  symmetry.  The  solutions 
presented  here  were  obtained  on  a  numerical  grid 
with  spatial  resolution  (in  the  flame  region)  better 
than  0.02  mm,  which  is  adequate  to  resolve  flames 
in  the  H2/air  system.  We  utilized  accurate  time  in¬ 
tegration  to  simulate  the  actual  procedure  employed 
in  the  experiment  to  produce  multiple  vigorously 
burning  flames. 


2-D  DNS  OF  OPPOSED-JET  FLAME 


803 


Fig.  1.  Isocontours  of  velocity  components,  tempera¬ 
ture,  and  mole  fraction  of  selected  species.  Left  column,  a 
contiguous  diffusion  flame;  right  column,  an  edge  flame. 
Re  =  508.8  in  both  cases;  also,  the  numerical  domain  is 
truncated  in  the  radial  direction.  Because  of  the  axial  sym¬ 
metry,  only  half  of  the  numerical  domain  is  shown  (the 
numerical  domain  on  the  right  is  the  mirror  image  of  that 
on  the  left). 

The  geometry  was  that  of  an  axially  symmetric  op¬ 
posed-jet  burner,  with  fuel  and  air  nozzles  of  diam¬ 
eter  D  =  0.27  cm  placed  at  a  distance  of  L  =  0.27 
cm  apart.  The  numerical  domain  was  large  (the  out¬ 
flow  boundary  was  placed  three  diameters  away 
from  the  nozzles)  to  accommodate  the  edge  flame, 
which  was  found  at  positions  away  from  the  axis  of 
symmetry.  The  geometry  of  the  numerical  domain 
was  simply  a  rectangle,  similar  to  that  reported  in 
Ref.  [24],  The  upper  and  lower  boundaries  of  the 
numerical  domain  (Fig.  1)  were  assumed  to  be  non¬ 
reactive  walls  of  constant  temperature  (fixed  at  300 
K).  The  flux  of  species  through  these  walls  was  zero, 
and  the  no-slip  boundaiy  condition  was  enforced. 
Zero-Neumann  boundary  condition  was  applied  at 
the  outflow  boundary.  We  are  aware  that  in  the  ex¬ 
periment  reported  in  Ref.  [1],  there  was  no  cold 
wall,  and  air  entrainment  would  further  dilute  the 
mixture  in  the  boundary.  We  chose  a  cold  wall  as  our 
numerical  boundary  in  order  to  obtain  a  well-defined 
numerical  system  and  to  avoid  the  need  to  model 
the  effects  of  air  entrainment.  However,  it  is  a  simple 
matter  to  add  two  plates  in  future  experiments  to 
match  the  numerical  domain  utilized  here.  The 
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Fig.  2.  Flame  structure  of  the  diffusion  flame  found 
along  the  axis  of  symmetry  at  Re  =  508.8. 

boundary  conditions  (at  nozzle  exits)  were  as  follows: 
“parabolic”  (and  thus  non-uniform)  velocity  profiles 
were  prescribed  at  the  exits  of  the  fuel  and  oxidizer 
nozzles  with  equal  peak  velocities  of  300  cm/ s  (cor¬ 
responding  to  Re  =  508.8  using  the  fluid  properties 
of  air  at  300  K).  Fuel  consisted  of  0.4  H2  and  0.6  N2 
(by  volume)  at  300  K,  while  the  oxidizer  was  air  at 
300  K.  These  physical  parameters  were  chosen  from 
one  set  of  the  experimental  results  reported  in  Ref. 
[1],  in  which  multiple  vigorously  burning  solutions 
and  hysteresis  were  observed.  In  subsequent  sec¬ 
tions,  all  results  are  non-dimensionalized  by  the  peak 
velocity  of  the  airstream  at  the  inflow  boundary  (U^r 
=  300  cm/s)  and  the  diameter  of  the  nozzles,  D. 
Time  is  non-dimensionalized  by  tc  =  D/U^. 

The  steady-state  two-dimensional  diffusion  flame 
from  which  we  begin  our  calculations  was  first  ob¬ 
tained  as  in  Ref.  [24] .  We  first  simulated  the  mixing 
of  fuel  and  oxidizer  in  the  axisymmetric  two-dimen¬ 
sional  opposed-jet  setup,  then  applied  an  ignition 
source  and  computed  the  entire  ignition  process  un¬ 
til  all  transients  died  out  and  the  steady  two-dimen¬ 
sional  diffusion  flame  (at  Re  —  508.8)  was  estab¬ 
lished. 


Results  and  Discussion 

We  focus  on  two  different  steady  solutions — a 
classical  diffusion  flame  and  an  edge  flame;  both 
were  obtained  under  identical  conditions.  The  disk¬ 
shaped  diffusion  flame  at  Re  =  508.8  is  depicted  in 
Fig.  1  (left  column).  Because  of  the  small  size  of  the 
nozzles,  their  close  placement  relative  to  one  an¬ 
other,  and  the  parabolic  outflow  velocity  boundary 
condition,  the  strained  diffusion  flame  considered 
here,  strictly  speaking,  cannot  be  accurately  de¬ 
scribed  by  the  commonly  used  one-dimensional 
model  [29,30],  The  flame  structure  found  along  the 
axis  of  symmetry  is  shown  in  Fig.  2.  This  diffusion 
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outflow 


Fig.  3.  Isocontours  of  the  heat  release  rate, 
(biHi)/ (/)Cp),  which  is  non-dimensionalized  by  tc  and  300 
K;  superimposed  on  this  plot  are  the  streamlines,  the  stoi¬ 
chiometric  line,  and  two  linear  coordinates  (lines  a  and  b), 
along  which  the  flame  structures  are  extracted. 
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Fig.  4.  Flame  structure  of  the  diffusion  flame  (tempera¬ 
ture  and  mole  fractions  of  H2  and  02)  found  along  the  front 
of  the  edge  flame  (along  line  a,  as  indicated  in  Fig.  3). 
Distance  is  non-dimensional. 


flame  was  subjected  to  a  near-extinction  scalar  dis¬ 
sipation  rate;  a  further  increase  of  the  Reynolds 
number  extinguished  the  flame  in  the  region  near 
the  axis  of  the  system.  The  transition  Reynolds  num¬ 
ber  agrees  with  what  was  reported  in  Ref.  [1]. 

Starting  with  the  steady  diffusion  flame  obtained 
at  Re  =  508.8,  the  Reynolds  number  was  increased 
to  700,  and  the  diffusion  flame  disk  was  first  extin¬ 
guished  in  the  vicinity  of  the  axis  of  symmetry.  An 
edge  flame  formed  and  propagated  radially  outward 
until  it  was  stabilized  at  a  certain  radial  position, 
forming  a  new  steady  flame.  The  stabilization  mech¬ 
anism  seen  here  is  aerodynamic  in  nature  due  to  the 
divergence  of  the  velocity  field  in  the  radial  direc¬ 
tion.  Starting  with  this  stable  solution  (an  edge 


flame)  at  Re  =  700,  the  Reynolds  number  was  de¬ 
creased  to  Re  —  508.8.  The  edge  flame  propagated 
back  toward  the  axis  (an  ignition  wave)  but  stabilized 
at  a  finite  distance  from  the  axis.  The  original  disk 
shape  diffusion  flame  did  not  recur  even  though  the 
scalar  dissipation  rate  had  been  decreased  to  its 
original  subcritical  level.  Isocontours  of  velocity 
components,  temperature,  and  the  mole  fractions  of 
selected  species  of  these  steady  flames  are  depicted 
in  Fig.  1.  The  left  column  corresponds  to  the  diffu¬ 
sion  flame  and  the  right  column  to  the  edge  flame. 
Note  that,  for  clarity,  the  numerical  domain  has  been 
truncated  in  the  radial  direction  in  the  figure. 

Of  interest  is  the  structure  of  this  edge  flame, 
which  propagates  against  the  flow  of  a  non-uniform 
mixture.  In  regions  ahead  of  the  edge  flame,  the 
mole  fraction  of  H2  varied  in  the  range  of  4%-17% 
(in  the  “preheat”  region  slightly  ahead  of  the  edge 
flame).  To  further  analyze  this  edge  flame,  we  plot¬ 
ted  in  Fig.  3  the  isocontour  of  the  heat  release  rate 
((ZiWjH^/ipCp))  non-dimensionalized  by  fc/300  K. 
We  also  plotted  the  streamlines  and  the  stoichio¬ 
metric  line  (Z  =  Zst).  Here,  Z  =  (yYHo  —  Yq2  + 
Yo2.air)/(yYH2,fuel  +  Yo,,air)>  where  y  denotes  the 
stoichiometric  ratio — that  is,  y  =  1/2W0,/ WH„  (Wj 
being  the  molecular  weight  of  species  i ).  The  stoi¬ 
chiometric  mixture  fraction  was  Zst  =  0.3892  for  the 
fuel  and  oxidizer  compositions  considered.  Several 
interesting  features  can  be  seen  in  these  figures: 

1 .  The  stoichiometric  mixture  fraction  line  does  not 
go  through  the  region  where  the  heat  release  rate 
is  at  its  maximum  value.  Instead,  it  passes  near 
the  cold  wall,  where  chemical  reaction  is 
quenched  and  the  edge  flame  only  has  a  lean  arm. 

2.  The  slanted  orientation  of  the  edge  flame  and  the 
leakage  of  hydrogen  (marked  by  a  circle  in  Fig. 
1)  are  seen  near  the  cold  wall  boundary.  Hydro¬ 
gen  diffuses  toward  the  high-temperature  zone 
and  is  fully  consumed  in  regions  downstream  of 
the  edge  flame.  However,  there  is  no  discernible 
diffusion  flame  structure  downstream  of  the  edge 
flame. 

3.  The  streamlines  diverge  as  they  approach  the 
flame  front  (consistent  with  the  results  of  Ref. 
[16]). 

To  further  illustrate  the  flame  structure,  we  extract 
cross  sections  (in  directions  “parallel”  and  “orthog¬ 
onal”  to  the  flame  front)  along  the  lines  shown  in 
Fig.  3.  The  origins  (black  dots)  and  directions  of  the 
one-dimensional  coordinates  on  which  the  data  are 
reported  in  Figs.  4  and  5  are  also  shown  in  Fig.  3. 
Fig.  4  is  taken  along  the  line  approximately  parallel 
to  the  edge  flame  in  the  preheat  zone  (line  a  in  Fig. 
3).  The  plots  in  Fig.  4  indicate  clearly  that  the  edge 
flame  sees  a  mixture  of  hydrogen  and  air  with  strong 
variations  in  concentration  along  the  direction  par¬ 
allel  to  the  flame  front.  A  rich  mixture  exists  in  the 
small  region  near  the  wall  boundary  on  the  fuel  side 


2-D  DNS  OF  OPPOSED-JET  FLAME 


805 


Fig.  5.  Flame  structure  of  the  diffusion  flame  (tempera¬ 
ture;  mole  fractions  of  H2,  OH,  and  02;  velocity  projected 
on  to  the  linear  coordinate,  and  non-dimensional  heat  re¬ 
lease  rate)  found  orthogonal  to  the  front  of  the  edge  flame 
(along  line  b,  as  indicated  in  Fig.  3).  Distance  is  non-di¬ 
mensional. 

(the  right  edge  in  Fig.  4).  However,  the  cold  wall 
there  quenches  the  flame. 

The  plots  of  temperature,  species  concentrations, 
and  the  velocity  component  orthogonal  (projected 
along  line  b  in  Fig.  3)  to  the  edge  flame  are  reported 
in  Fig.  5.  Note  that  line  b  in  Fig.  3  goes  through  the 
point  where  the  heat  release  rate  attains  its  maxi¬ 
mum  value,  and,  around  this  region,  the  plotted  pro¬ 
files  in  Fig.  5  resemble  those  of  a  premixed  flame. 
However,  far  upstream  and  downstream  of  the  re¬ 
action  zone,  the  scalar  profiles  show  significant  var¬ 
iations  reflecting  the  two-dimensional  nature  of  this 
edge  flame.  For  example,  it  can  be  seen  that,  down¬ 
stream  of  the  reaction  zone,  02  concentration  rises 
and  temperature  drops  due  to  transport  in  the  trans¬ 
verse  direction.  The  propagation  velocity  of  the  edge 
flame  (taken  to  be  the  velocity  of  the  flow  ahead  and 
in  a  direction  orthogonal  to  the  flame  front)  is  about 
105  cm/ s,  which  agrees  well  with  that  of  a  lean  1 12/ 
air  mixture  [31], 

The  edge  flame  is  responsible  for  the  hysteresis 
reported  in  Ref.  [1].  After  the  diffusion  flame  disk  is 
extinguished,  an  edge  flame  that  has  the  character¬ 
istics  of  a  premixed  flame  forms.  This  edge  flame 
stabilizes  itself  at  a  location  where  the  incoming  flow 
velocity  matches  the  local  propagation  velocity.  To 
recover  the  initial  contiguous  diffusion  flame  ( con¬ 
tiguous  refers  to  a  simply  connected  flame  as  op¬ 
posed  to  one  with  a  hole),  one  has  to  reduce  the  flow 
rate  such  that  the  local  propagation  velocity  of  this 
edge  flame  exceeds  the  local  flow  velocity  to  allow  it 
to  propagate  all  the  way  back  to  the  axis  of  symmetry. 
Reducing  the  flow  rate  to  produce  a  subcritical  strain 
rate  alone  is  not  sufficient  to  ensure  that  a  contigu¬ 
ous  diffusion  flame  can  be  recovered. 


Conclusions 

Using  detailed  chemistry  and  realistic  transport, 
we  performed  direct  numerical  simulations  to  ex¬ 
plain  the  transition  between  two  different  vigorously 
burning  flame  configurations  observed  in  a  recent 
experiment.  We  found  two  stable  solutions  that  ex¬ 
plain  the  multiple  solutions  described  briefly  in  Ref. 
[1].  Our  DNS  results  revealed  the  nature  of  the  sec¬ 
ond  solution:  it  is  an  edge  flame  that  propagates 
against  a  stream  of  non-uniform  mixture  in  the  op¬ 
posed-jet  configuration.  The  hysteresis  observed  in 
the  experiment  can  be  explained  by  the  propagation 
of  the  edge  flame.  This  hysteresis  is  different  from 
those  associated  with  the  classical  S-shaped  curve: 
for  the  same  value  of  the  Reynolds  number,  two  so¬ 
lutions  corresponding  to  vigorously  burning  flames 
coexist  with  the  pure  mixing  (or  weakly  burning 
flame)  solution.  R  may  also  be  important  in  improv¬ 
ing  our  understanding  of  the  local  reignition  in  a 
turbulent  flame  brush,  in  that  the  reduction  of  the 
local  strain  rate  to  the  subcritical  level  may  not  be 
sufficient  for  the  diffusion  flame  to  re-establish  due 
to  the  propagation  character  of  the  edge  flame.  The 
local  velocity  field  has  to  allow  the  propagation  of 
this  edge  flame  to  reignite  the  mixture  and  re-estab¬ 
lish  the  contiguous  diffusion  flame.  Future  work  will 
aim  at  utilizing  an  arc-length  continuation  method 
to  compute  the  solution  branches  of  this  transition 
phenomenon  and  to  construct  experiments  that 
match  the  numerical  system  more  closely.  Lastly,  the 
agreement  between  the  numerical  results  and  what 
was  observed  in  the  experiment  suggests  that  the 
traditional  opposed-jet  experimental  setup  can  be 
used  as  a  device  to  study  edge  flames.  One  clear 
advantage  of  this  setup  is  that  the  system  is  strictly 
two-dimensional,  making  it  possible  to  have  match¬ 
ing  detailed  numerical  simulations  to  complement 
the  experimental  endeavor. 
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COMMENTS 


Jaideep  Ray,  Sandia  National  Labs,  USA.  Did  you  ob¬ 
serve  the  tribrachial  flame  in  both  the  one-step  and  H2/02 
cases?  Also,  if  you  plotted  the  fuel  consumption  rate  of  the 
tribrachial  configuration,  was  it  symmetrical  about  the  stoi¬ 
chiometric  line? 

Author’s  Reply.  We  did  not  perform  calculations  of  H2/ 
02  flames  with  one-step  chemistry.  In  our  earlier  paper  [1], 
the  calculation  which  showed  tribrachial  flame  structure 
utilized  a  fictitious  fuel/oxidizer  system  with  one-step 
chemistry  and  simplified  transport.  The  transport  proper¬ 
ties  and  the  stoichiometry  (1:1)  there  does  not  correspond 
to  the  H2/02  system;  it  was  merely  a  generic  reactive  sys¬ 
tem  used  to  illustrate  our  numerical  approach  which  was 
new  at  the  time.  Nevertheless,  when  we  changed  the  Dam- 
kohler  number,  this  generic  system  (on  an  opposed-jet  con¬ 
figuration)  changed  from  the  traditional  diffusion  flame 


disk  to  a  tribrachial  flame  structure.  Because  of  the  sim¬ 
plification  adopted  in  this  earlier  work,  the  tribrachial  flame 
was  completely  symmetrical  on  both  the  fuel  and  oxidizer 
sides  (all  thermodynamic  properties  of  the  fuel  and  oxidizer 
are  identical  and  the  mass  fractions  of  both  reactants 
streams  are  equal  to  one  at  the  nozzle  exits).  We  would  like 
to  note  that  in  this  simple  calculation,  when  the  stoichi¬ 
ometry  was  changed  from  1  to  0.5,  the  tribrachial  flame 
becomes  slanted  in  much  the  same  way  as  the  H2/02  edge 
flame  shown  in  the  realistic  simulation  reported  in  the  cur¬ 
rent  endeavor. 
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